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ABSTRACT 

Detailed chemical abundances are presented for seven M31 outer halo globular clusters (with 
projected distances from M31 greater than 30 kpc), as derived from high resolution integrated 
light spectra taken with the Hobby Eberly Telescope. Five of these clusters were recently 
discovered in the Pan-Andromeda Archaeological Survey (PAndAS)—this paper presents 
the first determinations of integrated Fe, Na, Mg, Ca, Ti, Ni, Ba, and Eu abundances for 
these clusters. Four of the target clusters (PA06, PA53, PA54, and PA56) are metal-poor 
([Fe/H] < — 1.5), a-enhanced (though they are possibly less a-cnhanccd than Milky 

Way stars at the ler level), and show signs of star-to-star Na and Mg variations. The other 
three globular clusters (H10, H23, and PA17) are more metal rich, with metallicities ranging 
from [Fe/H] = —1.4 to —0.9. While H23 is chemically similar to Milky Way field stars. 
Milky Way globular clusters, and other M31 clusters, H10 and PA 17 have moderately low 
[Ca/Fe], compared to Milky Way field stars and clusters. Additionally, PA17’s high [Mg/Ca] 
and [Ba/Eu] ratios are distinct from Milky Way stars, and are in better agreement with the 
stars and clusters in the Farge Magellanic Cloud (EMC). None of the clusters studied here 
can be conclusively linked to any of the identified streams from PAndAS; however, based on 
their locations, kinematics, metallicities, and detailed abundances, the most metal-rich PAn¬ 
dAS clusters H23 and PA 17 may be associated with the progenitor of the Giant Stellar Stream, 
H10 may be associated with the SW Cloud, and PA53 and PA56 may be associated with the 
Eastern Cloud. 

Key words: galaxies: individual(M31) — galaxies: abundances — galaxies: star clusters: 
general — globular clusters: general — galaxies: evolution 


1 INTRODUCTION 

The Andromeda Galaxy (M31) is the nearest large neighbour to 
the Milky Way (MW). Its relative proximity offers a unique op¬ 
portunity to study the stellar populations of a nearby, large galaxy 
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without being hindered by interstellar extinction. Though the stars 
in M31 are considerably fainter than MW field stars, they can still 
be resolved for photometric and low resolution spectroscopic ob¬ 
servations while entire globular clusters (GCs) can be targeted for 
integrated light (IL) spectroscopy. Such observations contribute to 
studies of the formation and evolution of M31 and its satellite sys¬ 
tem. 

Deep surveys of M31’s o uter regions, such as the Isaac 
Newt on Telescope (INT) survey jlbata et al.ll200ll : iFerguson et al.l 
2002|) and the Pan-Andro meda Archaeological Survey (PAndAS; 
McConnachie et ahl 120091) have demonstrated that there is a sig¬ 
nificant amount of coherent substructure in the outer halo; these 
coherent streams indicate that M31’s halo is currently being as- 
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sembled by accretion of dwarf satellites. Ilbata et ail j2014h present 
maps of the outer halo, demonstrating that there are a multitude of 
streams, plumes, and clouds of metal-poor stars, along with a sig¬ 
nificant metal-rich population that is mostly located in the Giant 
Stellar Stream (GSS) south of M31. As much as 42% of the most 
metal-poor stars ([Fe/H] < —1.7) lie in coherent substructure— 
for metal-rich stars ([Fe/H] > —0.6) this percentage rises to 86%. 
The rest of the outer halo stars are distributed in a “smooth” com¬ 
ponent. 

The presence of metal-rich and metal-poor coherent streams 
implies t hat M31 is currently accreting multiple dwarf satellites 
(see, e.g. Jjohnston et alj2008l for accretion simulations that lead to 
coherent streams). The locations and metallicities of the outer halo 
field stars can be used to infer the nature of these dwarf galaxies. 


(1) The metal-rich GSS and the surrounding metal-rich features 
indicate that M31 has accreted at lea s t one fairly massive 
(~ LMC-mass) galaxy (j lbata et all200l|;|Ferguson e t al.l2002l ; 


Ilbata et al.l2007ljFardal et all2013l : Ibata et alj 20141) . Star for¬ 

mation histories indicate that most of the metal rich stars 
formed more than 5 Gyr ago, hinting that the prog enitor may 
have been an early-type dwarf ( [Bernard et al.ll2014f ). 

(2) Intermediate [Fe/H] stars are als o found in cohere n t substruc¬ 
tures, notably the SW Cloud ( Ilbata et al.1120141 ; iBate et all 
20141). which has its own GCs (PA7, PA8, and PA 14; 
Mackey et al.1 20131. |2014|) and possibly its own H I gas 


dLewis et al .1 120 1 3F Estimates based on metallicity and total 
brightness suggest that the SW Cloud progenitor was also a 
fairly massive dwarf gala xy (w ith a mas s comparable to the 
Fornax dwarf spheroidal; IBate et al.ll2014li . 

(3) The most metal-poor streams are gas-free (based on 21 cm 
observations; iLewis et alj[2013l ) and likely originated in lower 
mass dwarf spheroidal systems. 


The outer halo GCs also suggest that M31 has had a fairly 
active accretion history. Of the ~ 60 GCs discovered in PAn- 
dAS, many appear to lie along stellar streams. Mackey et al .1 j2010l) 
demonstrated that the positions of these GCs are correlated with the 
positions of the streams, with a < 1% chance that these GCs are lo¬ 
cated along the streams by chance. PA7 and PA8 are kinematically 
associated with each other 1 Macke^_etalJ2013b and with the SW 
Cloud itself dBate et alj2014 : lMackev et al.l2014l) . IVelianoski et al.1 
J2013L l2014h also show that many of the GCs which appear to 
lie along streams could be associated, based on their kinematics. 
IColucci et al.112014) have also demonstrated that at least one outer 
halo GC (G002) has likely been accreted from a dwarf satellite. 

Taken as a whole, the entire outer halo GC system is also dif¬ 
ferent from the MW’s GC population. For instance, iHuxor et al.l 
l i2014l) show that M31 has more luminous and faint GCs than the 
MW—they suggest that M3 l_could have acquired these extra GCs 
through accretion. lVelianoski et al.l ( 12014b) also demonstrate that the 
entire GC system has clear signs of rotation in both the GCs spa¬ 
tially associated with streams and those that are unassociated with 
streams. This rotation suggests either 1) that M31 experienced a 
merger with a galaxy large enough to bring in a substantial fraction 
of the outer halo GCs, or 2) that the parent dwarf galaxies were ac¬ 
creted from a preferred direction such that their angular momenta 
were correl ated, as is pres ently seen for a separate plane of dwarf 
satellites ( Ilbata et al.ll2013t) . 

Ultimately, the observations of outer halo stars and GCs indi¬ 
cate that M31 has recently experienced a merger with at least one 
massive dwarf and multiple lower mass dwarfs. Detailed chemi¬ 
cal abundances (from high resolution spectroscopy) are an effec¬ 


tive way to isolate and identif y coeval groups, a process known 
as “chemical tagging ” (e.g., iFreeman & Bland-Hawthorr] 120021 ; 
iMitschane et al.ll2014l) . With abundances of a, iron-peak, and neu¬ 
tron capture elements, metal rich ([Fe/H] > —1.5) stars and 
GCs from dwarf gala xies can be distinguished from those in mas¬ 
sive g alaxies (e.g. IVenn et ai1l2004l : iTolstov et al .1120091 ; iTing et ail 
120121) . Unlike GCs associated with the MW and its dwarf satellites, 
however, distant extragalactic targets must be observed through 


be determined from lower resolution IL spectra (e.g. 

Schiavon et al. 

2002 

; Lee & Worthev|2005:; Puzia et al.12006. 200$ 

Caldwell et al. 

2011 

), the detailed chemical abundances necessary for chemical 


tagging analyses require high resolution IL spectr oscopy (see, e.g., 


|McWj]liam&Bernstem[ 201)4 SaktMetldJ 20J3|)^As in the MW 


dCohenll2004l ; ISbordone et al ■1200.4 ISakari et alJl201 lh . if the M31 
outer halo GCs originated in lower mass dwarf galaxies then they 
could have distinct abundance patterns from MW field stars, de¬ 
pending on metallicity. 

High resoluti on IL spectroscopy has been t e sted extensively 
on Galactic GCs dMcWilliam & Bemsteinll2008l : Ic amer otj 120091 ; 
ISakari et al .1120 1 3I . I20I4I) . demonstrating that 


(1) IL abundances trace the individual stellar abundances and rep¬ 
resent cluster averages when abundances do not vary between 
stars 

(2) Abundance determinations from high resolution IL spec¬ 
troscopy are more precise than from lower resolution studies 

(3) Certain abundance ratios are very stable to uncertainties in the 
underlying stellar populations (e.g. interloping field stars, un¬ 
certainties in cluster age, microturbulence relations, etc.), such 
as [Ca/Fe] 

(4) Partially resolved photometry down to the horizontal branch 
(HB) reduces systematic uncertainties by constraining [Fe/H] 
and HB morphology. 


Thus, a high resolution IL spectroscopic analysis of M31 GCs 
can provide precise and accurate abundances that will be suit¬ 
able for chemical tagging. High resolution IL spectroscopic tech- 
nique s have been applied to n earby GCs in the Local Group (in 
M31. IColucci et al.|2009l.l20H the LMC. IColucci et al.l201 lal and 
IColucci et alj 20121 ; and in dwarf galaxies, Coluccietak[ 201 lbli 


and in the large elliptical galaxy NGC 5128 ( Colucc i et al 


2013). 


At the distance of M31, the size and brightness of its GCs make 
them ideal targets for high resolution IL spectroscopy. 

This paper presents a high resolution IL chemical tagging 
analysis of seven outer halo M31 GCs, five of which were first 
observed in PAndAS. This is the first chemical analysis for these 
five PAndAS clusters, at any resolution. Section[2]presents the data 
and the analysis methods, while Section [3] describes the process 
for generating synthetic stellar populations for each cluster. The IL 
abundances are given in Section[4] and the implications for the for¬ 
mation of M31’s outer halo are then discussed in Section]?] Finally, 
the results are summarized in Section [6] 


2 OBSERVATIONS AND ANALYSIS METHODS 
2.1 Target Selection 

The primary goal of this high resolution spectroscopic investiga¬ 
tion of M31 GCs is to study the nature of the stars and clusters in 
M31 ’s outer halo. The targets were therefore restricted to GCs with 
projected distances from the centre of M31 that are 7? pro j > 30 
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kpc; priority was given to the GCs with the largest projected radii. 
Observational constraints required that the targets were sufficiently 
bright to observe in a reasonable amount of time and to avoid 
stochastic sampling issues—as a result, the target GCs are on_the 
brigh t end of the M31 outer halo GC distribution tsee lHuxor et ali 
l20l4 Seven outer halo GCs were targeted; their properties are 
summarized in TableQ] Tw o of these targ et clusters (H10 and H23) 
were discovered by iHuxor et al.l J2008h: the other fiv e were dis¬ 
covered in the PAndAS programme iHuxor et alj[20l4) . The loca¬ 
tions of these GCs on a density map of metal-poor outer halo stars 
(from M. Irwin, private communication) is shown in Figure|T] Note 
that none of the targets fall precisely onto any of the major stel¬ 
lar streams or overdensities in the outer halo. Nonetheless, PA56 
sits on a faint north-western extension of the Eastern Cloud, and is 
likely associated with this feature (Mackey et al. 2015, in prep.). 
While H23 has been tentatively linke d to S tream D based on its rel¬ 
ative proximity in projection dMackev et alj20 1 7)1). its r adia l velo c- 
ity r enders a genuine as sociation unlikely dVelianoski et al]|2014l) . 

ISakari et al.l d2014h present detailed tests of systematic abun¬ 
dance errors that occur because of uncertainties in the underlying 
stellar populations. These tests include: 

• Uncertainties in isochrone parameters such as age and [Fe/H] 

• Properties of evolved horizontal branch (HB) and asymptotic gi¬ 
ant branch (AGB) stars, particularly the HB morphology 

• Assumptions about the stellar populations (e.g. microturbulence 
relations, initial mass functions, mass segregation) 

• Unusual stars like interloping field stars, long period variables, 
etc. 


spectra, and variance weighting was used during aperture extraction 
to remove cosmic rays. Sky spectra (from the separate sky fibres) 
were replaced with continuum fits with the emission lines added 
back in (according to the sky line identifications from the UVES 
Quality Control sky spectrum websit^3) and were then subtracted 
from the object spectra. Telluric standards were also observed for 
removal of atmospheric absorption features. 

In order to avoid removing broad, blended features, the tar¬ 
get spectral orders were normalized with continu um fits to a n 
extremely metal-poor star (CS29502-092; see ISakari et al Jl20 13h . 
Low-order polynomial fits were then necessary to fully normal¬ 
ize the target spectra. The individual observations were cross- 
correlated with a high r esolution , high S/N Arcturus spectrum 
(from the Arcturus Atlafl; II in kkft20031) to determine radial veloc¬ 
ities. Heliocentric velocities were determined for each individual 
observation; the final, averaged heliocentric velocities are shown in 
Table [2] The quoted errors represent the dispersion between obser¬ 
vations. 

The individual, rest-frame spectra were combined with av¬ 
erage sigma-clipping rejection routines to remove any remaining 
cosmic rays; each spectrum was weighted by flux during the com¬ 
bination. The combined spectra were once again cross-correlated 
with the Arcturus template spectrum to determine cluster disper¬ 
sions (see ISakari et alj|2013l for a description of how this is done). 
These velocity dispersions are also shown in Table[2] Examples of 
the final, combined spectra are shown in Figure[2] 

2.3 Line List and EW Measurements 


These errors can be mitigated if GC properties are constrained 
with photometry or lower resolution analyses over a broader wave¬ 
length range. Observing priority was therefore given to those tar¬ 
gets with partially resolved photometry. Six of the seven target GCs 
have Hubble Space Telescope ( HST) photometry down to the HB 
dMackev et alJl2007l .Mackev et al. 2015, in prep.). This photome¬ 
try is used to place constraints on the best fitting isochrones (see 
Section[3]i. 


2.2 Observations and Data Reduction 

The targets were ob s erved with the Hobby -Eberly Telescope (HET; 
iRamsev et alJI 1 9981 : IShetrone et alj|200% at McDonald Observa¬ 
tory in Fort Davis, TX in 2011 and early 2012. The High Res¬ 
olution Spectrograph (HRS; fTuIIll 19981) was utilized with the 3" 
fibre and a slit width of 1", yielding an instrumental spectral res¬ 
olution of R — 30, 000. With the 600 gr/mrn cross disperser 

set to a central wavelength of 6302.9 A, wavelength coverages of 
~ 5320 — 6290 A and ~ 6360 — 7340 A were achieved in the 
blue and the red, respectively. The 3" fibre provided coverage of the 
clusters past their half-light radii; the additional sky fibres (located 
10" from the central object fibre) provided simultaneous observa¬ 
tions for sky subtraction. Exposure times were calculated to obtain 
a total S/N = 80 (per resolution element), although not all targets 
received sufficient time to meet this goal. The details of the obser¬ 
vations are shown in Table[2] 

Data reduction was performed in in the Image Reduction and 
Analysis Facility program (IRAF)Q As in lSakari et al.l d2013l ). bias 
removal was not performed to avoid adding noise to the individual 


1 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 


For these chemical tagging analyses, s pectral l ines o f Fe, Na,_Mg, 
Ca, Ti, Ni, Ba, and Eu are utilized. ISakari et~a D d20li l2014h 
demonstrated that for Galactic GCs reliable Fe, Ca, Ti, Ni, and 
Ba abundances can be determined from equivalent widths (EWs), 
while Na, Mg, and Eu abundances can be determined with spec¬ 
trum syntheses. Note that in this analysis the Ba lines are synthe¬ 
sized instead of using EWs, due to the lower S/N ratios of the target 
spectra. 

EWs for Fe, C a, Ti, and Ni lines were m e asured with the pro- 
gram DAOSPECfl l lStetson & Pancin(1l2008l) . ISakari et al .1 d201 3h 
verified that DAOSPEC measurements of IL spectral lines compare 
well with measurements from other methods for GCs with a range 
of velocity dispersions. Given the low S/N of some of the PAndAS 
targets, many of the DAOSPEC EWs were verified or refined by 
hand. IL spectral lines have contributions from stars with a vari¬ 
ety of line strengths—in particular, IL lines with EWs >100 mA 
can have contributions from bright red giant stars with EWs ~ 200 
mA. These str ong lines are dif ficult to model (see the discussions 
in Paper I and McWill iam et alj|~1995bl) and are typically removed 
from individual stellar analyses. For this reason, as a conservative 
estimate, all lines stronger than ~ 110 mA were removed from this 
analysis unless absolutely necessary. 

Table[3]shows the spectral lines that were used to derive abun¬ 
dances in the analysis, their atomic data, and the measured IL EWs 
for each cluster. The line list is based on the standard IL lines lists 


omy, Inc., under cooperative agreement with the National Science Founda¬ 
tion. 

2 http://www.eso.org/observing/dfo/quality/UVES/pipeline/sky_spe< 

3 ftp://ftp.noao.edu/catalogs/arcturusatlas/ 

4 DAOSPEC has been written by P.B. Stetson for the Dominion Astro- 
physical Observatory of the Herzberg Institute of Astrophysics, National 
Research Council, Canada. 
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Table 1. Properties of the target PAndAS clusters. 


Cluster 

RA (J2000) 

Dec (J2000) 

Lint 

no 

^proj 

(kpc) 

Photometric b 
[Fe/H] 

Spatial Association 
with Stream? 

References 

H10 

00 h 35 m 59.7 s 

-t-35°41'03".6 

15.7 

78.5 

-1.84 

N 

1.2 

H23 

00 h 54 m 25.0 s 

+39°42'55".5 

16.8 

37.1 

-1.54 

Stream D? 

1.2 

PA06 

00 h 06 m 12.0 s 

+41°41'21".0 

16.5 

93.7 

MP 

N 

2, 3,4 

PA17 

00 h 26 m 52.2 s 

+38°44'58".l 

16.3 

53.9 

?C 

N 

2, 3 

PA53 

01 h 17 m 58.4 s 

+39° 14'53". 2 

15.4 

95.9 

MP 

N 

2, 3,4 

PA54 

01 h 18 m 00.1 s 

+39° 16'59".9 

15.9 

95.8 

MP 

N 

2, 3,4 

PA56 

01 h 23 m 03.5 s 

+41°55'11".0 

16.8 

103.3 

MP 

N 

2, 3,4 


References: 1 = lMackey et alj i2Q07h , 2 = IVelianoski et alj fc014l) , 3 = lHuxor et alj J2014I) , 4 = Mackey et al., in prep.. 
a Projected distances are from the centre of M31. 

b Photometric metallicity estimates are used to constrain the parameters of the input isochrones; see Section [3 [Fe/H] estimates are based on Galactic GC 
fiducial fits. MP indicates a metal-poor cluster ([Fe/H] < —1.5); fiducial fits will be given in Mackey et al., in prep.. 
c PA 17 does not have an HST CMD, and there is therefore no a priori information about its metallicity or HB morphology. 
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Figure 1. Locations of the target PAndAS clusters on a metal-poor density map of the full PAndAS footprint. The locations of stellar streams and satellite 
galaxies are obvious as dark regions. The red stars show the metal-rich ([Fe/H] > —1.5; see Section[4} PAndAS clusters observed in this paper, while the 
blue stars show the metal-poor GCs. Note that PA53 and PA54 are very close together in projection and share a single point in the plot. 
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of iMcWilliam & Bernstein! i2008h and lColucci et al.l ll2Q09li . with 
supplements from the R ed Giant Br a nch (R GB) stellar line lists of 
I Sakari et al .1 J20 1 lh and IVenn et ail d2012l) . Lines that were mea¬ 
sured with spectrum syntheses are indicated. The full line lists to 
synthesize all lines in_a 10 A w i ndow (for spectrum syntheses) 
are the ones used in ISakari et al .1 d2013l) : they consist of the EW 


line list, with supplements from the Vienna Atomic Line Databas^EI 
(VALD: lKupka et alj2QQ0l l. the Kurucz database ^and the National 
Institute of Standards and Technology (NISTJJ database. Molec¬ 
ular features from the Kurucz database were included when in- 


5 http://www.astro.uu.se/~vald/php/vald.php 

6 http://kurucz.harvard.edu/linelists.html 

7 http://www.nist.gov/index.html 
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Table 2. PAndAS cluster observations. 


Cluster 

Observation 

Dates 

Exposure 
Time (s) 

S/N a 
(5500 A) 

S/N a 
(7000 A) 

r 5 

'ils 

(r h ) 

^helio,obs 

(km/s) 

^helio,lit 

(km/s) 

^obs 

(km/s) 

H10 

2011 Jan 2, 10. 11.22. 23,28, 30 

19180 

82 

140 

2.3 

—351.9±1.5 

-3523=9 

6.63=0.4 

H23 

2011 Jul 5, 7, 10. 11, Aug 2, 4, Sep 23 

16050 

65 

82 

4.1 

-373.33=0.1 

-3773=11 

6.23=0.4 

PA06 

2011 Sep 29, Oct 1, 2, 4, 6. 17, 18, 19, 21, 23 

29278 

65 

105 

2.6 

-341.43=0.7 

-3273=15 

5.63=0.4 

PA17 

2012 Jan 16, 18,21 

8100 

41 

50 

2.6 

-260.03=1.0 

-2793=15 

6.13=0.5 

PA53 

2011 Aug 1, Sep 20, 26 

8100 

115 

148 

2.7 

-270.83=0.9 

-2533=10 

12.03=0.4 

PA54 

2011 Aug 25. Sep 24. 25,27 

10800 

82 

130 

2.2 

-344.93=0.8 

-3363=8 

7.53=0.4 

PA56 

2011 Oct 7, 17, 18, 20, 21, 23, 24, 30, 

30623 

65 

82 

2.4 

-241.43=1.7 

-2393=8 

6.43=0.4 


Nov 19, Dec 30, 
2012 Feb 10, 11 


References: Literature radial velocities are from lVelianoski et alj 1 201 4) 

a S/N ratios are per resolution element, and assume th at there are 2.7 pixel s pe r resolution e l ement for HRS. 

b The coverage radii are based on the half-light radii in lTanvir~et •aim and lHuxor et al m2oh — all GCs are covered past their half-light radii. 


dicated in the Arcturus Atlas (for CH, CN, and MgH). Isotopic 
and hyperfi n e stru cture (HFS) components for the Ba II lines are 
from McWilliaml I 1 9981) while the Eu II com ponents are from 
lLawler et~ 1 2001 all and lLawler et al.l d2001bt) . The HFS compo¬ 
nents are included in the syntheses; while they do not significantly 
affect the strengths of the lines, they can affect the line profiles. 


2.4 Solar Abundances 

ISakari et aD d20 1 3l. l2014h prese nted measurements of the solar 
spectrum from the lKuruc j 12005!) solar flux atlas@ The [Fe/H] and 
[X/Fe] ratios in this paper were calculated line by line using the 
solar abundances derived from those EWs and spect rum s_yntheses. 
When the solar lines were stronger than 150 mA. lAsplund et al.l 
J2009I) solar abundances were adopted for those lines. 


are then assigned corre sponding Kurucz model atmosphere^ 
dCastelli & Kuruc j|2004l) . where the grid values are interpolated to 
a specific box’s T e ff and log g. When the [Ca/Fe] ratio is supersolar, 
o-enhanced model atmospheres are used. Microturbulent velocities 
are determined via an empirical relation b etween £ (in km s -1 ) and 
log q that fits the Sun and Arcturus (see IMcWilliam & Bernstein! 
120081 ) . 


The abundances of the Fe I lines (from the EWs in Table 
[3} are used to co nstrain the best-fitting isochrone’s single age and 
metallicity. A s in McWilliam & Bemsteinl I 2 OO 8 I ) and lColucci et al.l 
J2009I . |2011 al . |20 1 3^ . the adopted isochrone is the one that leads to 
the flattest trends in Fe I abundance with wavelength, reduced EW 
(REW)Fl and excitation potential (EP); th ese crit eria are adopted 
from analyses of individual stars (see IMcWilliam & Bernsteinl 

lHH)). 


3 UNDERLYING STELLAR POPULATIONS 


The first step in an IL analysis is to determine the atmospheric 
parameters of the stars in the underlying populations. The cen¬ 
tral regions of the PAndAS GCs cannot be resolved with HST, 
and high quality CMDs cannot be obtained for the regions in¬ 
cluded in the IL spectra—isochrones must therefore be used to 
model the stellar population s. Synthetic populations are generated 
with the BaSTI isochrones dPietrinfemi et alj|2004 ICordier et al.l 
1 20071) with exte nded AGBs an d mass loss parameters of 7 / = 
0.2, assuming a Kroupal (|2002|) IMF and the total cluster mag¬ 
nitudes from IHuxor et al.l | |2014|) . The BaSTI isochrones are uti- 
lized becaus e they extend all the way through the AGB phase— 
ISakari et al .1 d2014l) show that differences between isochrone sets 
are insignificant for the abundances derived in this paper. The 
synthetic Hertzsprung Russell Diagram (HRD) populations are 
binned into boxes, each with 3.5% of the total flux. The boxes 


8 http://kurucz.harvard.edu/sun.html 


3.1 Clusters with HST Photometry 

The stars in the outer regions of six of the targets clusters have 
been resolved with HST, enablingconstraints to be placed on the 
GC [Fe/H] and HB morphologvl 11 ! This photometry is not used 
to identify the isochrone that best fits the resolved CMD, but is 
only used to constrain the possible isochrones and HB morphology. 
These constraints help to drastically reduce the systematic errors in 
the IL abundances (see lSakari et al.|[2014l) . For the blue HB clus¬ 
ters (PA06, PA53, PA54, and PA56; see Figure d, synthetic BaSTI 
HB0 are used instead of the default HBs. ISakari et aD | |2014|) 
show that as long as the HBs are modelled approximately correctly 
the integrated abundances will not be systematically affected. 

The parameters of the spectroscopically-determined 


9 http://kurucz.harvard.edu/grids.html 

10 REW = EW/A 

11 Rough constraints can also be placed on cluster age, specifically 
whether the cluster is older than ~ 4 Gyr; additionally, the presence of 
RR Lyrae stars would constrain the age to > 10 Gyr. 

12 http://basti.oa-teramo.inaf.it/BASTI/WEB_TOOLS/HB_SYNT/index. 1 
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Figure 2. IL spectra of the PAndAS GCs, arranged by metallicity. Spectral lines used in typical RGB stellar analyses are identified. 


isochrones are listed in Table [4] and the isochrones them¬ 
selves are shown on top of the observed CMDs in Figure [3] (for 
the six clusters with HST photometry). Again, these isochrones 
were not derived based on fits to the observed CMDs, but were 
determined by minimizing trends in the Fe I line abundances. The 
fact that the isochrones agree well with the CMDs confirms the 
validity of these high resolution IL techniques. 

The slopes of the trends in Fe I abundances are given in Ta¬ 


ble [5] Again, solutions were selected by simultaneously minimiz¬ 
ing trends in Fe I abundance with wavelength, REW, and ER All 
REW and EP trends for the PAndAS clusters are flat within their 
lcr errors. The line-to-line dispersion can still be quite large, which 
leads to some uncertainty in isochrone age (~ 5 Gyr), though 

this h as a minimal effect on the integrated abundances l lSakari et al.1 
12014 . Also note that the BaSTI isochrones have large spacings in 
metallicity (~ 0.2 — 0.3 dex, depending on [Fe/H]), preventing 
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Table 3. The Line List. a 


A 

(A) 

Element 

E.P. 

(eV) 

log gf 

Equivalent width (mA) 

H10 H23 PA06 PA17 

PA53 

PA54 

PA56 

5324.191 

Fe I 

3.21 

-0.103 

- 

- 

80.0 

83.1 

79.4 

80.0 

5339.937 

Fe I 

3.27 

-0.720 

81.0 

- 

- 

- 

- 

54.8 

5367.476 

Fe I 

4.42 

0.443 

77.9 

93.9 

46.0 

68.0 

50.0 

53.2 

5369.974 

Fe I 

4.37 

0.536 

- 

- 

50.0 

- 

- 

59.4 

5371.501 

Fe I 

0.96 

-1.644 

- 

- 

126.0 

- 

- 

- 


Notes: Equivalent widths were measured in DAOSPEC; all strong lines were checked and refined in splot. IL lines stronger than 110 mA were not included 
in the analysis. 

“Table [3] is published in its entirety in the electronic edition of Monthly Notices of the Royal Astronomical Society. A portion is shown here for guidance 
regarding its form and content. 


Table 4. Parameters of the “best-fitting” HRDs and synthetic and observed 
integrated colours. 



Isochrone 

[Fe/H] 

Age 

(Gyr) 

Synthetic 
(V - J)o 

Observed 0 

(V - 1)o 

A(V-/) 0 

H10 

-1.31 

12 

0.94 

0.95 

-0.01 

H23 

-1.31 

9 

0.98 

1.01 

-0.03 

PA06 

-1.84 

12 

0.91 

0.87 

+0.04 

PA17 6 

-1.01 

12 

1.05 

1.14 

-0.09 

PA53 

-1.62 

12 

0.89 

0.85 

+0.04 

PA54 

-1.84 

13 

0.89 

0.84 

+0.05 

PA56 

-1.62 

12 

0.88 

0.89 

-0.01 


“ Uncertainties in the observed colours from lHuxor et aft 1 20 111) are likely 
to be ~ 0.05 dex. 

b PA17 does not yet have a high-quality / 1ST CMD; therefore, there are 
few a priori constraints on the underlying stellar population. 

exact metallicities from being adopted for the isochrones. Compar¬ 
isons with the integrated colours are shown in Table [4] The differ¬ 
ences for the six partially resolved GCs are all < 0.05 magnitudes 
(whic h is t h e appr oximate uncertainty in the observed colours from 
iHuxor et al.ll2014l) : this indicates that the populations are at least 
reasonably well modelled. H10 and H23’s slightly redder synthetic 
colours may be due to the fact that the default HBs are slightly too 
red (see Figures [3al and [3bt. while the inconsistencies in the other 
GCs might be due to incorrect ages or slight discrepancies in the 
RGBs, HBs, or AGBs. 


3.2 PA17: The Cluster without HST Photometry 

Since PA 17 does not yet have a high quality HST CMD there is 
little a priori information available to constrain its age, [Fe/H], or 
HB morphology. As a result the isochrone default HBs were used. 
The flattest slopes occur for a relatively metal-rich isochrone, with 
[Fe/H] ~ —1, and an age of 12 Gyr. Given PA17’s [Fe/H] and red 
integrated colour ((V — I) o = 1.14, which is redder than H10 and 
H23; IHuxor et aPl2014h , the default BaSTI HB (which is red) is 
likely to be sufficient, and there abundance offsets as a result of the 
uncertain HB morphology should be negligible. With the default 
HBs, the Fe I abundances imply a best-fitting age and [Fe/H] of 
12 Gyr and -1.01, respectively, for PA17. 


3.3 Comparisons with Literature Values 


HlO’s isochrone age and int egrated [Fe/Hl agree well with the 
high resolution IL values from lColucci et all I I2014I) . which are de¬ 
termined from spectra with a slightly different wavelength range. 
However, the metallicities and ages for H10 and H23 from this 
work do not agree with the previous photometric results from 
iMackev et al.l d2007l). From analys e s of p artially resolved CMDs 
of the upper RGBs, IMackev et al.l d2007h find [Fe/H] values for 
H10 and H23 that are lower than the integrated [Fe I/H] ratios 
presented here. The fiducial fits in Mackey et al. are performed 
with an optimization routine that finds the best combination of 
Galactic GC fiducial (at a given [Fe/H]), distance modulus, and 
reddening that fit the observed RGB. This should provide a good 
estimate of cluster [Fe/H], since the slope o f the RGB is pri¬ 
marily sensitive to metallicity (e.g. ISaraiedinill 19941) , though this 
[Fe/H] is dependent on the adopted distance modulus and redden¬ 
ing. For instance. [Mackey et all d2007h find that HlO’s RGB falls 
between the fiducials of M92 ([Fe/H] = -2.14) and NGC 6752 
([Fe/H] = —1.54), implying that H10 likely has a metallicity of 
[Fe/H] ~ —1.84. However, Figure l4al demonstrates that HlO’s 
RGB slope can be well represented by M3’s RGB, indicating a 
higher metallicity of [Fe/H] ss —1.5 (which agrees with the spec¬ 
troscopic values from this work). 

A similar technique can be applied to all the PAndAS GCs. 
An empirical calibration between Galactic GC RGB slope (in the 
F606W and F814W HST filters) is provided in Appendix [A] The 
RGB slopes and metallicities of the PAndAS GCs are shown in 
Figure[4b]on top of this empirical Galactic GC relation. The cluster 
HBs were first aligned with 47 Tuc and/or M3 to fix the distance 
modulus and reddening. The approximate (V — I) o colours at two 
fixed absolute magnitudes (My = 0 and My = — 2) were then 
determined, and RGB slopes were calculated. The approximate un¬ 
certainties in RGB slope are estimated from the uncertainties in 
RGB colour and distance modulus—note that PA56 has a large 
uncertainty because its paucity of bright HB stars makes its dis¬ 
tance modulus more uncertain. The quoted errors i n [F e/Hl a re 0.1 
dex, to reflect possible systematic errors JSakari et al.ll2014l ). The 
PAndAS GCs agree with the Galactic GC relation within the la 
errors; however, all the PAndAS GCs have slightly steeper slopes 
than their Galactic counterparts, hinting that the clusters may be a- 
deficient compared to the calibrating MW GCs (see the discussion 
in Appendix[A}- Regardless, this general agreement shows that the 
spectroscopic metallicities derived in this paper are consistent with 
the observed CMDs. 
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Table 5. Trends in Fe I abundance for the PAndAS Clusters 


Cluster 

Wavelength 

Slope 

10 -5 dex/A 

REW 

Slope 

(dex) 

EP 

Slope 

(dex/eV) 

H10 

2.7±4.0 

-0.030i0.097 

—0.011±0.019 

H23 

- 1.5±2.9 

—0.016±0.067 

-0.0020i0.015 

PA06 

- 1.8±2.3 

-0.0057i0.058 

0.0041i0.011 

PA 17“ 

3.3±7.5 

0.035±0.19 

0.015i0.029 

PA53 

3.7±5.7 

0.019±0.18 

0.019i0.025 

PA54 

9.3±3.7 

—0.00018±0.099 

—0.020i0.020 

PA56 

6.8±5.6 

0.033i0.13 

0.0081i0.026 


a PA17 does not yet have a high-quality HST CMD; therefore, there are few a priori constraints on the underlying stellar population. 


4 ABUNDANCES 

Abund ances were determined with ILABUNDS (first introduced in 
McWilliam & Berns teinl2008l) . an IL modification of the 1997 ver¬ 
sion of the Local Thermo dyna mical Equilibrium (LTE) line anal¬ 
ysis code MOOG ( Snedeniri973h . The final IL abundance^ for all 
elements are shown in Table[6] R andom errors for EW -based abun¬ 
dance s are calculated as in lShetrone et al.l J2003I) and I Sakari et al .1 
J201 lh . where the largest of three errors was adopted: the line-to- 
line scatter for that element, the error due to EW uncertainties, or 
the iron line-to-line scatter. The largest of these three uncertainties 
is adopted as the final random abundance error for that element. 
Random errors for abundances derived via spectrum syntheses are 
calculated based on u ncerta inties in continuum placement and line 
profile fitting (see lSakari et al.ll20l3l) . 

Systematic errors are adopted from lSakari et al .1 J20l4) . based 
on GC [Fe/H] and HB morphology. The optimal abundance ratios 
for chemical tagging are those which 

• are least sensitive to uncertainties in the underlying stellar popu¬ 
lations 

• are most distinct between massive and dwarf galaxies. 

Sakari et al. determined that the most stable and useful abundance 
ratios were [Ca I/Fe I], [Ni I/Fe I], and [Ba II/Eu II]. The first 
two [X/Fe] ratios are typically stable (within ~ 0.1 dex) to most 
systematic uncertainties in the underlying populations. Though in¬ 
dividually [Ba II/Fe II] and [Eu II/Fe n] are highly sensitive to 
the underlying stellar populations, the Ba II and Eu II offsets are 
often similar, and the [Ba/Eu] ratio is much more stable to pop¬ 
ulation uncertainties (in a partially resolved system the uncertain¬ 
ties can still be as high as ~ 0.1 dex). These ratios are also use¬ 
ful for comparisons with dwarf galaxy stars. The [Ca/Fe] ratio 
can serve as an [a/Fe] indicator, since metal rich dwarf galaxy 
stars have lower [Ca/Fe] at a given [Fe/H] (see, e.g.. iTolstov et al J 

13 The standard notation is used, i.e. 

[X/H] = log e(X)—log e(X) 0 = (l2+log - (l2+log 
where TVx is the column density of any element X. 


120091) . The [Ni/Fe] ratio has been observed to be low in metal- 
rich field stars an d clusters associated with the Sagit t arius (Sgr) 
dwarf spheroidal dCohenl 120041 : ISbordone et all 120051120071) and 
anomalous MW field stars that are suspected to h ave been ac¬ 
creted from dwarf galaxies dNissen & Schusterll2010t) . Finally, the 
[Ba/Eu] ratio is o bserve d to be high in dwarf galaxy stars (see, 
e.g„ iTolstov et alj|2009 i). Other suitable, high precision ratios in¬ 
clude [Eu/Ca] and [Mg/Ca]. Although the other element ratios 
([Na I/Fe I], [Mg I/Fe I], [Ti I/Fe I], and [Ti II/Fe II]) are less stable 
to uncertainties in the underlying populations, their agreement (or 
disagreement) with other abundance ratios may prove useful. 

Table [7] presents estimates of the maximum systematic er¬ 
rors in the integrated abundances of the PAndAS clu sters. All er¬ 
rors estimates are based on the tables presented in 1 Sakari et al J 
d2014l) . with the exception of Na and Mg, which wer e calculated 
in the same way. The errors are selected from I Sakari et aid d2014l) 
to match the metallicity and HB morphology of the cluster, as well 
as the analysis technique. The GCs H10, H23, PA06, PA53, PA54, 
and PA56 are partially resolved, which reduces the systematic er¬ 
rors considerably. At [Fe/H] ~ — 1.4, H10 is closest in [Fe/H] to 
M3, M13, and NGC 7006; its systematic uncertainties are therefore 
averages from the three clusters (though note that since all HlO's 
HB stars are red, the hot star uncertainties are not included). H23 
has a metallicity in between 47 Tuc and M3/M13/NGC 7006, and 
its quoted uncertainties are therefore an average of the 47 Tuc val¬ 
ues and the mean M3/M13/NGC 7006 values. PA53, PA54, and 
PA56 are similar in [Fe/H] and HB morphology to M3, M13, and 
NGC 7006, and their averaged abundance offsets are therefore uti¬ 
lized. M15’s systematic uncertainties are adopted for PA06. For 
PA 17, the unresolved cluster, the systematic uncertainties of 47 Tuc 
were adopted (since PA17’s [Fe/H] ~ —0.9 is closest to 47 Tuc’s 
[Fe/H] ~ —0.8; see lSakari et all2013i) . assuming that the cluster 
is modelled with isochrones and that there are no constraints on the 
underlying population. 

In the following subsections, each element type is discussed 
separately. The [Fe/H] and [X/Fe] ratios of the PAndAS GCs are 
plotted along with field stars in the MW and its dwarf satellites. 
There are no detailed abundances for M31 field stars, though there 
are medium resolution a-abundances for four outer halo M31 stars 
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(a) H10 (b) H23 (c) PA06 





(d) PA53 (e) PA54 (f) PA56 

Figure 3. HST CMDs of the six partially resolved PAndAS GCs. The spectroscopically determined best-fitting isochrones are shown (in red) on the CMDs to 
illustrate how the synthetic HRDs are populated. For the blue HB clusters PA06, PA53, PA54, and PA56, synthetic HBs (shown in blue) were used instead of 
the default HBs. 


l lVargas et al.ll2014h : however, M31 field star abundance patterns 
are likely similar to the MW, especially for metal poor stars. The 
PAndAS clusters are also compared to integrated and averaged 
abundances from GCs in the MW, the inner halo of M31, and MW 
satellite galaxies. 

4.1 Iron 

The [Fe/H] abundances were determined with EWs, as described in 
Section 12731 Because 1) there are fewer F e II l ines and 2) Fe II is 
more susceptible to systematic effects dSakari et al.l2013j) . [Fe I/H] 
is adopted as the representative [Fe/Fl], despite potential non-LTE 
(NLTE) effects. Note that for all clusters Fe I and Fe II are in rea¬ 
sonably good agreement. Three clusters, H10, H23, and PA17, are 
more metal-rich than [Fe/H] = — 1.5, while the other four are 

metal-poor. The most metal-poor clusters all have fewer detectable 
lines (~ 30 Fe I lines and 1-2 Fe II lines) than the more metal-rich 
GCs, leading to larger random errors in [Fe/H]. 

Most GC systems in dwarf galaxies tend to be metal-poor, 
with more massive dwarf galaxies possessing more metal-rich GCs 


(e.g. IPeng et all 120061) . The comparatively high metallicities of 
H10, H23, and PA17 therefore suggest that they formed in a galaxy 
more massive than the Fornax (For) dwarf spheroidal. It is clear 
from Tables[]]and[6]that the metal-rich PAndAS GCs have smaller 
projected distances from the centre of M31 than the more metal- 
poor GCs. Based on a l ow reso lution survey of inner halo, bulge, 
and disk M31 GCs. ICaldwell et ail j201 lh found evidence for an 
abundance gradient in the inner regions of M31—this gradient 
seems to flatten to [Fe/H] ~ —1.8 at distances > 2 kpc. H23 and 
PA 17 are therefore more metal-rich than expected given their large 
distance from the centre of M31; however, metal-rich field stars 
have been identified in streams in M31’s outer halo dlbata et alj 
» 

4.2 o-elements 

The a-elements (Mg, Ca, and Ti) form primarily through cap¬ 
tures of 4 He nuclei during hydrostatic burning in massive stars, 
while Fe forms both during core collapse supernovae of massive 
stars and during the detonation of white dwarfs whose progeni- 
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Figure 4. Left: The H10 HST CMD, plotted on M3’s V, I CMD from the ACS Globular Cluster Survey jSaraiedini et alj|2007l ; I Anderson et alj|2008b . The 
distance modulus and reddening have been adjusted to match the horizontal branches of the two clusters. The slopes of the GCs agree well, supporting the 
higher spectroscopic metallicity from this study. Right: RGB slopes (determined from the HST photometry) versus spectroscopic metallicity. The solid line 
shows the calibration with Galactic GCs (presented in Appendix 0, while the red stars show the slopes of the PAndAS clusters. The RGB slopes predict 
metallicities that agree well with the integrated [Fe/H] values. The blue dot-dashed line shows the offset that would occur at a fixed [Fe/H] if clusters have 
[a/Fe] = 0 (see Appendix lAl. 


Table 6. Integrated abundances of the PAndAS clusters, with random errors and the number of lines for each element. 



PA17 

H23 

H10 

PA53 

PA56 

PA54 

PA06 

[Fe I/H] 

-0.93 ±0.03 

-1.12 ±0.02 

-1.36 ±0.02 

-1.64 ±0.03 

-1.73 ± 0.03 

-1.84 ±0.02 

-2.06 ±0.02 

N 

31 

65 

45 

29 

40 

32 

30 

[Fe II/H] 

-0.96 ±0.15 

-1.21 ±0.08 

-1.30 ±0.04 

-1.68 ±0.10 

-1.71 ± 0.10 

-1.75 ±0.08 

-2.09 ±0.10 

N 

4 

5 

5 

1 

1 

2 

1 

[Na I/Fe I] 

0.60 ±0.15 

-0.01 ±0.23 

0.35 ±0.25 

0.41 ±0.30 

0.43 ± 0.20 

0.78 ± 0.30 

0.71 ± 0.25 

N 

2 

2 

1 

1 

2 

2 

2 

[Mg I/Fe H 

0.80 ±0.15 

0.30 ±0.20 

0.43 ±0.15 

0.20 ±0.20 

0.54 ±0.15 

0.32 ± 0.18 

0.14 ±0.20 

N 

1 

2 

2 

1 

2 

2 

1 

[Ca I/Fe I] 

0.04 ±0.07 

0.41 ±0.04 

0.25 ±0.04 

0.19 ±0.03 

0.24 ±0.05 

0.28 ± 0.04 

0.46 ±0.07 

N 

4 

6 

10 

10 

8 

9 

9 

[Ti I/Fe I] 

0.0 ±0.10 

0.56 ±0.05 

0.39 ±0.10 

0.16 ±0.10 

0.23 ±0.10 

- 

- 

N 

1 

3 

1 

1 

1 

0 

0 

[Ti II/Fe II] 

0.17 ±0.10 

0.39 ±0.10 

0.40 ± 0.03 

0.17 ±0.06 

_ 

0.49 ± 0.10 

0.31 ± 0.10 

N 

1 

1 

2 

2 

0 

1 

2 

[Ni I/Fe I] 

-0.08 ±0.09 

-0.07 ±0.04 

0.01 ±0.04 

-0.11 ±0.18 

-0.10 ±0.16 

-0.05 ±0.06 

-0.02 ±0.06 

N 

3 

5 

7 

2 

2 

3 

2 

[Ba II/Fe II] 

0.21 ± 0.25 

-0.03 ±0.11 

-0.05 ±0.11 

0.03 ±0.22 

-0.39 ± 0.15 

-0.10 ±0.20 

-0.17 ±0.20 

N 

2 

2 

2 

2 

1 

1 

1 

[Eu II/Fe II] 

0.36 ±0.35 

0.53 ±0.25 

0.75 ±0.20 

0.68 ±0.25 

0.96 ±0.15 

< 0.50 

< 1.04 

N 

1 

1 

1 

1 

1 

1 

1 


Notes: [Fe/H] and [X/Fe] ratios are calculated differently, line by line, relative to the solar abundances for that line (see I Sakari et ali2013l and Sakari et alj 
120141 for the EWs used to derive the solar abundances). 
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Table 7. Systematic error estimates. 


A [Fe/H] 


A[X/Fe] 


A[X/Y] 



Fe I 

Fe II 

Nal 

Mg I 

Cal 

Til 

Ti II 

Ni I 

Ball 

Eu II 

[Mg/Ca] 

[Na/Mg] 

[Ba/Eu] 

[Eu/Ca] 

PA17 

0.19 

0.21 

0.21 

0.27 

0.13 

0.16 

0.17 

0.07 

0.23 

0.13 

0.16 

0.15 

0.13 

0.11 

H23 

0.13 

0.19 

0.13 

0.14 

0.08 

0.12 

0.10 

0.07 

0.18 

0.10 

0.08 

0.11 

0.09 

0.09 

H10 

0.09 

0.15 

0.11 

0.07 

0.05 

0.08 

0.08 

0.07 

0.09 

0.08 

0.04 

0.10 

0.07 

0.10 

PA53, PA54, PA56 

0.11 

0.16 

0.13 

0.08 

0.05 

0.08 

0.11 

0.08 

0.10 

0.09 

0.05 

0.11 

0.11 

0.10 

PA06 

0.20 

0.13 

0.10 

0.08 

0.07 

0.11 

0.09 

0.07 

0.16 

0.10 

0.07 

0.10 

0.12 

0.12 


Notes: Total errors are conservatively estimated by adding the individual errors from lsakari et alJ ( 2014 ) in quadrature, based on the appropriate [Fe/H] and 
HB morphology. 


tors were lower mass stars. Comparisons between the abundances 
of cr-elements and Fe vs. [Fe/H] provide an indication of how the 
chemical contributions from different types of stars have changed 
in a specific environment. Because the trends in [a/Fe] with [Fe/H] 
differ between massive and dwarf galaxies, the abundance ratios 
can be used to link stars and/or GCs to their birth environments. 


4.2.1 Magnesium 


Magnesium abundances for the PAndAS GCs are presented in Ta¬ 
ble [6] For most clusters [Mg/Fe] is determined from the 5528 and 
5711 A lines—however, PA17’s [Mg/Fe] is only determined from 
the 5711 A line because its 5528 A line is too strong; similarly, 
the 5711 A line is too weak in the most metal-poor cluster, PA06, 
and in the GC with the highest velocity dispersion, PA53. Com¬ 
parisons with MW, M31, and dwarf galaxy field stars and clus¬ 
ters are shown in Figure [5] Field stars in the most massive, nearby, 
well-studied dwarfs that host GCs are also shown, including the 
LMC and the Sgr and For dwarf spheroidals; lower mass galax¬ 
ies such as Sculptor and Carina are not included since they do 
not have GCs. The average and/or integrated abundances of MW, 
M31, LMC, Sgr, and For GCs are shown separ ately from the field 
stars. The unusual MW GCs Palomar 1 (Pal 1; Sakari et ~al]|201 lh 
and Rupr echt 106 (Rup 106; Villanova et al.ll20 131) and the M31 
GC G002 IColucci et alJ d2014 ) are given different symbols, since 
their chemical abundances indicate that they may have originated 
in dwarf galaxies (though their host galaxies have not yet been iden¬ 
tified). 

Note that star-to-star Mg variations have been observed 
in ma ssive, metal- poor Galactic clusters ([Fe/H] < —1.2; 

ICarretta et al.1 l2009al) ; in those clu sters, Mg variatio n s can 
be as large as 0. 5 — 1 dex (e.g. ISneden et alJ 1 19971 1 20041 
ICohen & Melendezl 1 20051) . The most metal-poor PAndAS GC, 
PA06, has a low [Mg/F el ratio for its metallicity, similar to M15 
(see ISakari et al.l 120131) ; this suggests that strong star-to-star Mg 
variations may exist in PA06. The integrated [Mg/Fe] ratios in the 
other metal-poor GCs roughly agree with the G alactic GCs. The 
integrated Mg abundances of the IColucci et al.l i2014l) M31 GCs 
are systematically lower than the PAndAS GCs at the metal poor 
end—however, this may be due to the masses of the observed clus¬ 
ters. Many of the Colucci et al. targets have larger velocity disper¬ 
sions and higher total luminosities, indicating that they are more 
massive than these PAndAS GCs. If this is the case, the GCs in 
the Colucci et al. sample may harbour stro nger Mg variations (e.g. 
Carrett aet al.l2009al ).IColucci et al.l d2014l) do find a correlation be¬ 
tween integrated [Mg/Fe] and total My for the metal poor GCs, 
such that the more massive GCs have lower integrated [Mg/Fe]— 


all five of the most metal poor ([Fe/H] < —1.2) PAndAS GCs are 
consistent with the Colucci et al. relation for all GCs. 

Because of their high metallicities the two most metal-rich 
PAndAS GCs, H23 and PA 17, are not likely to host significant star- 
to-star variations in Mg, though they may host multiple popula¬ 
tions (see Section l4.3.2b . The [Mg/Fe] ratios of H10 and H23 agree 
with Galactic field stars . Note that HlO’s [Mg/Fe] is higher than the 
value presented by IColucci et al.l d2014l ). which was derived from 
slightly different lines; considering random and systematic uncer¬ 
tainties, however, the two are consistent. PA 17 has a [Mg/Fe] ra¬ 
tio of 0.8 dex, which is higher than most of the field stars and all 
of the GCs. The syntheses of PA17’s 5711 A line are shown in 
Figure [6]— this line has a similar strength in all 3 exposures, sug¬ 
gesting that the line strength is not due to an improperly removed 
cosmic ray. Note that there are no other Mg I features in the ob¬ 
served spectrum (other than the strong 5528 A line) to verify this 
high abundance. With its high [Mg/Fe], PA17 falls at the upper end 
of the Galactic and LMC field stars. This high [Mg/Fe] ratio sug¬ 
gests a greater contribution from massive star ejecta, either from 
a top-heavy IMF or from inhomogeneous mixing and stochastic 
sampling of supernova ejecta. Alternatively, high Mg could indi¬ 
cate contributi ons fro m a rapidly rotating massive star that went 
supernova fe.g. lTakahashi et al.ll2014l) . 


4.2.2 Calcium and Titanium 

Calcium and titanium mainly form explosively in core collapse su¬ 
pernovae and generally behave like the other a-ele ments, tho ugh 
both elements also form in Type la supernovae dNomoto et al .1 
11984 iMaeda et alj|2010h . The integrated [Ca I/Fe I] abundances 
are determined based on EWs of < 10 spectral lines. All Ca I 
line EWs were verified by hand, and any uncertain lines were re¬ 
moved from the analysis. Spectrum syntheses were also performed 
on several Ca I lines to verify the EW-based abundances. The Ti 
abundances are based on fewer lines (1-3), but both [ Ti I/Fe I] and 
[Ti II/Fe II] generally agree with [Ca/Fel. ISakari et all d2014l) de¬ 
termined that the [Ti/Fe] ratios are more sensitive to systematic un¬ 
certainties in the underlying stellar populations; thus, only [Ca/Fe] 
is compared with other targets. 

The PAndAS cluster [Ca/Fe] ratios are compared to MW 
and dwarf galaxy field stars in Figure [7a] and to MW, M31, and 
dwarf galaxy GCs in Figure [7b] Most of the PAndAS GCs are Ca- 
enhanced, though PA17’s [Ca/Fe] is roughly solar. All the PAn¬ 
dAS clusters look very similar to the LMC and Sgr GCs; most 
of the For clusters are more metal-po or tha n the PAndAS GCs. 
HlO’s [Ca/Fe] agrees well with lColucci et alJd2014l) and is slightly 
lower than the average MW and M31 clusters. PA17’s low [Ca/Fe] 
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Figur e 5. Top: Comparisons of PAndAS clusters (red stars) to MW field stars (grey; from the sources in I Venn et al.|2004i with supplements fromlReddy et al 
2006) and dwarf galaxy field stars. LMC stars ar e shown in blue JPompeia et alJ200 a). Sgr in magenta JSbordo ne et alJ2007l:l M onaco et alJ20071 : ICarretta et al 
201C ; IChou etlikkoiol : I Me William et al.ll2013l) , and For in green iTafelmever et alj|2010l ; IShetrone et al.ll2003t ; iLetarte et alj|2010l) . The error bars show the 
random errors. The dashed horizontal line shows the solar value. Random errors are shown as thick error bars, while systematic and random errors (added 
in quadrature) are shown with the thinner error bars. Bottom: Comparisons of PAndAS clusters (red stars) to MW field stars (grey), MW clusters (black 
circles), other M31 clusters (maroon), and dwarf galaxy GCs. IL abundanc es are shown as lar ger symbols, while averaged G C abundances from individual 
stars are shown as smaller sym bols. IL abundances of MW clusters are from l Sakari et alJ J2014I) , while cluste r averages are from l Pritzl et alJ i2005l) . The M31 
cluster IL abundance s are froml Colucci et alJ j2009tl2014l). LMC cluster s are in blue; IL abundances are from lColucci et al]fe012t), w h ile the individua l stella r 
abu ndances are fromljohnson et al] 120061) and Mucciarel li et alJ J2008l) . Average abundances of S gr clusters are fromlCohenl j2004l)JSbordone et alJ J 20051) , 
and lMottini et alJ J 20081). Averages of the individual stars in For clusters are in green, and are f rormLeta rte et al.| {20060 . The unusual MW clusters Palomar 1 
and Rup 106 (from lSakari et alJ[2oTll and lVillanova et aljf2013l) and the M31 GC G002 (from lColucci et alj|2014l7 are shown with crosses—these GCs may 
have been accreted from dwarf galaxies, though they have no obvious host galaxies. 


could be understood several ways. As a a-element, low [Ca/Fe] 
could be indicative of the slow s tar formation r atios or top-light 
IMFs seen in dwarf galaxies (e.g. iTins lev 1979; Me William et al .1 
|201 3l) — however, PA17’s high [Mg/Fe] ratio suggests that this in¬ 
terpretation is not so simple. Regardless of the cause, PA17’s Ca 
abundance looks most similar to Pal 1, Ter 7, Pal 12, and the LMC 
field stars and clusters. 


4.2.3 Mg vi. Ca 

If Mg and Ca have the same nucleosynthetic site then the [Mg/Ca] 
ratio should not change with [Fe/H]. In Galactic field stars this is 
more or less what is observed: [Mg/Ca] is roughly solar except at 


the lowest metallicities where the dispersion is very high. The dis¬ 
persion at low metallic ities is t ypically ascribed to inhomogeneous 
mixing (see, e.g., I Venn et all 120121) . However, metal-rich dwarf 
galaxy stars appear to have different [Mg/Ca] ratios from MW stars. 
In fact, the nucleosynthetic sites of Mg and Ca are not the same: 
Mg is m ainly produced during the hydrostatic burning of massive 
stars (e.g. lWooslev & Weaveil 1 9951) whil e Ca is formed explosively 
in both Type II and la supemovae (e.g. IWooslev & Weavetll 1 9951 . 
iNomoto et al.ll 1984L iMaeda et al.l l2010h . Figure [8] shows that the 
metal-rich field stars in the LMC and For have higher [Mg/Ca] 
ratios than MW stars at the same metallicity—this discrepancy is 
ty pic all y ascribed to the different nucleosynthetic sites of Mg and 
Ca (e.g. lShetrone et al]l2003l : IVenn et al]|2004l : ILetarte et alj|20ld : 
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Figure 6. Syntheses of the 5711 A Mg I line in PA17. The red line shows the best fit, while the green and blue lines are the dzlcr fits. 


IVenn et aill20 1 2l) . However, stars in Sgr hav e low [Mg/Ca] ratio s 
(except possibly f or star s in the Sgr streams: [Monaco et alJI 20071) ; 
[McWilliam et all 120131) attribute these low ratios to a top-light 
IMF, where the lack of massive stars means that less Mg is pro¬ 
duced than Ca. Most of the PAndAS GCs overlap with Galactic 
and dwarf galaxy stars; however, PA17’s [Mg/Ca] ratio is higher 
than the Galactic stars, in agreement with the LMC stars. 

Again, star-to-star Mg variations could affect the integrated 
[Mg/Ca] ratios of the metal-poor GCs, complicating comparisons 
with field stars. However, if the integrated [Mg/Fe] ratios are af¬ 
fected by multiple populations, one might expect the integrated 
[Mg/Fe] and [Mg/Ca] ratios to be lower than the “primordial” val¬ 
ues (because Mg is expected to decrease in the second popula¬ 
tion). Only PA06 has a lower [Mg/Ca] than the majority of the 
Galactic stars, in agreement with the integrated [Mg/Ca] of Ml5, 
which has known, s trong star-to-star Mg variations that af fect the 
integrated [Mg/Fe] ISneden et al.lll997l : lsakari et al.ll2013l) . Thus, 
PA06 seems to be similar to Ml5, and may host a significant sec¬ 
ond, Mg-deficient population. The other metal-poor GCs agree well 
with the Galactic stars and GCs, with the exception of PA56, which 
has a slightly high [Mg/Ca]. Note that most of the PAndAS GCs 
have higher [Mg/Ca] ratios than the M3 1 GC s—again, this is likely 
because the GCs in the lColucci et ali l20l4) sample are more mas¬ 
sive than the PAndAS GCs, and therefore may have stronger star- 
to-star Mg variations. 

The GCs more metal rich than [Fe/H] ~ —1.2 are not likely 
to host significant Mg variations. H10 and H23 have roughly so¬ 
lar [Mg/Ca] rati os, i n agreement jwith the Galactic field stars and 
clusters (though fColucci et ailhoidl find a lower [Mg/Ca] for H10 
because of their lower [Mg/Fe]). PA 17 is higher than the Galac¬ 
tic field stars and clusters by ~ 2 o, in agreement with the LMC 
and Fornax field stars. This high [Mg/Ca] is driven by the high 
[Mg/Fe] (determined from spectrum syntheses of the 5711 A line) 
and the low [Ca/Fe] (from EWs of four clean lines). Again, the 
high Mg suggests that PA 17 is enriched with the products of the 
most massive stars, which may be due to IMF effects, stochas¬ 
tic samplin g of supernova ejecta, or rotating supernova progeni¬ 
tors jTakahashi et alj2014l also note that rotating massive stars will 


not produce as much Ca, compared to non-rotating stars). Despite 
its agreement with the LMC field stars, PA17 is distinct from the 
LMC GCs, because the [Mg/Ca] ratios of the metal rich LMC GCs 
are not high like the LMC field stars. For the GC with the lowest 
[Mg/Ca] ratio (NGC 1718). [Colucci et all d2012h suggest an inho¬ 
mogeneous mixing scenario, where NGC 1718 formed out of mate¬ 
rial with no contributions from the highest mass stars. However, the 
three metal rich, [Mg/Ca]-poor LMC GCs have intermediate ages 
(1-3 Gyr), and it may not be entirely appropriate to compare these 
intermediate-age GCs with the older PAndAS GCs. Regardless of 
the disagreement with the LMC GCs, PA17’s high [Mg/Ca] ratio is 
most similar to the LMC and Fornax field stars. 

4.3 Sodium 

Sodium abundances are derived from spectrum syntheses of the 
6154 and 6160 A lines, which are all weaker than 100 mA in the 
PAndAS IL spectra. Sodium lines can be particularly sensitive to 
NLTE effects; however, the 6154/6160 A lines are not expected to 
have significant NLTE corrections (< 0.2 d e x) at subsolar metallic- 
ities and for weak line strengths fe.g. iLind et a 1 .1 [2oTT[> No cor¬ 
rections were applied to the [Na/Fe] abundances, which are shown 
in Table [6] 

4.3.1 [Na/Fe] 

Figure [9] shows comparisons between the PAndAS GCs and the 
field stars and GCs associated with the various galaxies. Several 
of the PAndAS GCs have [Na/Fe] ratios that are higher than MW 
and dwarf galaxy field stars. PA06, PA 17, and PA54 all have high 
[Na/Fe] ratios (> 0.6 dex) which places them above the Galac¬ 
tic and dwarf galaxy stars, while the H10, PA53 and PA56 ra¬ 
tios are mildly high (~ 0.4 dex). These high values agree well 
with the IL values of the Galactic and old LMC GCs and with 

14 Though note that some of the HRD boxes (see Section[3j can have EWs 
~ 130 mA in the most metal-rich clusters. 
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Figure 7. Top: Comparisons of PAndAS clusters (red stars) to MW field stars and dwarf galaxy field stars. Points and references are as in Figure[5a] Bottom: 
Comparisons of PAndAS clusters (red stars) to MW field stars and GCs from the MW, M31, and various dwarf galaxies. Points and references are as in Figure 
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many of the averaged Galactic GCs. In Galactic GCs the [Na/Fe] 
ratio is affected by star-to-star abundance variations within the 
clusters, whic h are often observed as Na/O anticorrelations (see 
ICarretta et alj|2009ah . These abundance variations are typically in¬ 
terpreted as arising from two chemically distinct stellar popula¬ 
tions within the GCs, which may indicate two se parate generations 
of stars (see the review by iGratton et alj|2012h . Several Galactic 
GCs have more centrally concentr ated Na- enh anced, O-deficient 
populations (e.g. M 13 and 47 Tuc: Ijohnson & Pilachowskill20 1 2| ; 
ICordero et al,l 120141). IL spectra of the cent ral regions (such as 
the Galactic GCs in lSakari et al]|20T3l . l2014t) will likely be dom¬ 
inated by the second generation (Na-enhanced, O-deficient) pop¬ 
ulations. It is therefore probable that (like the Galactic GCs) the 
PAndAS GCs are Na-enhanced because of the presence of multi¬ 
ple populations. It is possible that H23’s solar [Na/Fe] is due to 
the larger coverage of the IL spectrum (see Table |2), which may 
include light from the less centrally concentrated first generation 
sta rs; alternativ el y, H23 may have more “primordial” stars. Again, 
the lColucci et alj|2014l) [Na/Fe] ratio for H10 is lower than the one 
derived here—however, the ratios are derived from only two lines 


in both cases and are consistent within random and systematic er¬ 
rors. 

Thus, the [Na/Fe] ratios indicate that these PAndA S GCs are 
likelyto be “classical” GCs under the definition of ICarretta et al] 
(2009a . i.e. the GC stars exhibit the Na/O anticorrelation). Note 
that in this paper only Na abundances are derived for the PAndAS 
GCs, though O abundances would also be affected by multiple pop¬ 
ulations. However, the O lines are fairly weak in this wavelength 
range, and are not easily detectable in such low S/N spectra. Lines 
in other regions (e.g. the near infrared) would provide more robust 
integrated O abundances. 


4.3.2 [Nct/Mg] 

Na and Mg are expected to form primarily in massive stars (e.g. 
IWooslev & Weavenll995l : iTimmes et alJI 19951) : in the absence of 
multiple population effects, the abundances of Na and Mg are there¬ 
fore expected to correlate. The MW field stars with [Fe/H] > 
—2 sh ow a clear trend of increasing [Na/Mg] with [Fe/H]— 
iGehren et al] ]200d ) interpret this as a signature of increasing Na 
yields with metallicity. However, the metal poor (—2 < [Fe/H] < 
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Figure 8. Top: Comparisons of PAndAS clusters (red stars) to MW field stars and dwarf galaxy field stars. Points and references are as in Figure[5a] Bottom: 
Comparisons of PAndAS clusters (red stars) to MW field stars and GCs from the MW, M31, and various dwarf galaxies. Points and references are as in Figure 

[5b] 


— 1.2) Galactic and PAndAS GCs show the opposite trend: in PA06 
and M15 [Na/Mg] is very high, though [Na/Mg] decreases to sub¬ 
solar values by [Fe/H] ~ —1.5. The GCs therefore behave in a 
very different way from the field stars. Again, this is likely to be 
a result of multiple populations, combined with the metallicity de¬ 
pendency of the Mg/Al anticorrelation (e.g. lCarretta et al]l2009al ). 
For these lower mass GCs, at low metallicities the star-to-star Mg 
variations are likely to be strong in these (fairly bright) GCs, driv¬ 
ing the integrated [Mg/Fe] down and [Na/Mg] up. With increasing 
[Fe/H], the integrated [Mg/Fe] returns to its “primordial” value, and 
the [Na/Mg] ratio stabilizes. The [Na/Mg] abundances of the PAn¬ 
dAS GCs therefore also suggest that most of the GCs host multiple 
popul ations. Note th at the behaviour of the massive M31 clusters 
(from lColucci et alj|2014l) is similar at the metal poor end, though 
their [Na/Mg] ratios are higher at the metal rich end—again, this 
may be because the massive, metal rich GCs have stronger Mg vari¬ 
ations. 

PA17 remains an unusual object, because its similarly high 
[Na/Fe] and [Mg/Fe] ratios place its [Na/Mg] in agreement with the 
MW field stars, despite its disagreement in [Mg/Fe], If PA17’s high 
[Mg/Fe] is caused by increased supernovae ejecta from massive 


and/or rapidly rotating stars, then [Na/Fe] would also be enhanced; 
thus, PA17 may have an additional integrated Na enhancement that 
is caused by contributions from massive stars. 


4.4 Nickel 

Nickel is expec ted to fo rm in core collapse and Type la supernovae, 
along with Fe l lTimmes etal .11 19951) . The [Ni/Fe] ratios were de¬ 
termined with EWs, even though there are few (~ 1 — 7) Ni I 
lines—however, these lines were carefully checked to ensure that 
they were clean (i.e. uncontaminated by cosmic rays, sky lines, or 
noise) and were properly measured. The Ni abundances are given 
in Table [6] and are compared to field stars and other GCs in Figure 
EH The PAndAS clusters all have [Ni/Fe] ratios that are consistent 
with MW and dwarf galaxy stars—however, H23 and PA 17 both 
have slightly low Ni abundances, in better agreement with dwarf 
galaxy stars and GCs. In particular, PA17's [Ni I/Fe I] ratios agree 
well with the Pal 12 and Ter 7 ratios. 
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Figure 9. Top: Comparisons of PAndAS clusters (red stars) to MW field stars and dwarf galaxy field stars. Points and references are as in Figure[5a] Bottom: 
Comparisons of PAndAS clusters (red stars) to MW field stars and GCs from the MW, M31, and various dwarf galaxies. Points and references are as in Figure 
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4.5 Neutron Capture Elements: Ba and Eu 

The Ba and Eu abundances are determined through spectrum syn¬ 
theses. Only two Ba n lines were considered: the 5853 and 6141 A 
lines—the 6496 A line was removed because of p ossible NLTE 
effect s which lead to abundance overestimates dMashonkina et al] 
l200 0!). No molecular lines are included in these syntheses, since 
none are identified in those regions in the Arcturus Atlas. The Eu II 
abundances are determined through syntheses of the 6645 A lin e 
and include CH and CN molecular lines (see lSakari et al.ll2013l ): 
these syntheses are shown in Figure IT2l HFS and isotopic compo¬ 
nents were included for the Ba II and Eu II lines, as described in 
Section [2.3l 

I Sakari et al .112014l) demonstrated that the [Ba/Eu] ratio is of¬ 
ten more stable than the individual [Ba/Fe] and [Eu/Fe] ratios to 
uncertainties in the underlying population. Comparisons of [Ba/Eu] 
ratios are shown in Figure 1131 and illustrate that metal-rich dwarf 
galaxy stars have higher [Ba/Eu] ratios than MW field stars. PA06 
and PA54 only have upper limits on the Eu abundance, which trans¬ 
lates into a low er l imit in [B a/Eul. The estimated r-process-only ra¬ 
tio (from lBurris et al.l2000l) is also indicated in Figurefl3l The PAn¬ 
dAS GCs agree with the MW field stars and GCs, with the excep¬ 


tion of PA17 and PA56. PA56 has a low [Ba/Eu], much like M15. 
Both Ba and Eu have been known to vary within some of the most 
massive metal poor G Cf^l (such as M15; see, e.g.. ISneden et all 
I1997L lRoedererfl201 lh : neutron star mergers have been suggested 
as t h e origin of these star-to- star heavy element variations (e.g. 
iTsuiimoto & Shigevamal 120141) . PA56’s IL abundances may also 
be affected by the presence of heavy element dispersions with the 
GCs. Table [6] shows that PA56 has a moderate Ba abundance and 
a high Eu abundance, which leads to a [Ba/Eu] that is below the r- 
process only estimate; this indicates that PA56 may host stars with 
a significant heavy element dispersion. 

PA17’s high [Ba/Eu] is in agreement with the Galactic GCs 
within its 1 <t errors, though it agrees best with the dwarf galaxy 
stars and GCs, Pal 1, and Pal 12. The high [Ba/Eu] ratios in dwarf 
galaxies are typically interpreted as an excess of s-process over r- 
process elements. Its moderately high [Ba/Eu] ratio indicates that 
PA 17 has received chemical contributions from AGB stars. 

[Eu/q] is another popular chemical tagging indicator since 

15 Though note that not all massive GCs show heavy element variations, 
e.g. M92 lCohej201 ill . 
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Figur e 10. Top: Comparisons of PAndAS clusters (red stars) to MW field stars (grey; from the sources in Venn e t~alj2004l . with supplements from Re ddy ct alj 
1200611 and dwarf galaxy field stars. Points and references are as in Figure[5a] Bottom: Comparisons of PAndAS clusters (red stars) to MW field stars (grey), 
MW clusters (black circles), inner halo M31 clusters (maroon), and dwarf galaxy GCs. Points and references are as in Figure Bbl The dashed line shows the 
clear linear relationship for the metal poor GCs ([Fe/H] < —1.2). 


dwarf galaxy stars and GCs have higher [Eu/cr] ratios than MW 
stars and GCs at a given [Fe/H]—this h as been inter preted as a 
sign of an addition al r-process site ( e.g. Letarte et alj|20ld ) 


top-light IMF (e.g. iMcWilliam et al.1120131) . Figures 1741 and I14bl 
show the [Eu/Ca] ratios (which serve as [Eu/a] indicators) in differ¬ 
ent environments. Again, only upper limits are available for PA06 
and PA54. In general the clusters are in agreement with the MW 
and dwarf galaxy field stars and GCs, with the exception of PA56, 
which has a high [Eu/Ca], similar to M15—again, this is likely a 
signature of star-to-star Eu variations. 


5 DISCUSSION 

The detailed abundances of stars are dictated by the chemical com¬ 
position of the interstellar medium in their host galaxy. Chemi¬ 
cal comparisons between the PAndAS GCs and the field stars and 
GCs associated with other galaxies therefore provide clues about 
the nature of a GC’s birth environment, such as whether a cluster 
formed in a massive galaxy or a dwarf galaxy. Chemical compar¬ 
isons can also determine if the GCs are chemically similar to par¬ 


ticular galaxies, streams, or other GCs. If these GCs originated in 
dwarf galaxies, the abundances in this paper can therefore be used 
to infer the nature of the dwarfs that are currently being accreted 
into M31’s outer halo. 


5.1 Individual Clusters: A Summary 

5.1.1 PAndAS Clusters with [Fe/H] > —1.5 

The more metal-rich clusters all have higher metallicities than ex¬ 
pected given their large projected distances from the centre of M31. 

However, each GC has a unique chemical signature. 

PA17: Without HST photometry, PA17’s abundances are more 
uncertain than the other GCs; however, the high [Na/F el ratio 
sugges ts that PA 17 is a “classical” GC (under the ICarretta et al.1 
l2009al definition) with signs of a Na spread (and therefore 
probably an O spread as well). PA 17 is the most metal-rich of 
the target clusters, at [Fe/H] ~ —0.9. Despite the lack of HST 
photometry, systematic errors in [Fe/H] are likely to be < 0.2 
dex. PA17’s metallicity alone indicates that it formed in a fairly 
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Figure 11. Top: Comparisons of Ni abundances in the PAndAS clusters (red stars) to those in MW field stars (grey), dwarf galaxy field stars (coloured points), 
and the IL abundances of the Galactic clusters (black circles). Points and references are as in Figure [5] Bottom: Comparisons of [Ni/Fe] ratios in PAndAS 
clusters to MW field stars and clusters, Pal 1, and dwarf galaxy clusters. Points and references are as in Figure l5bl 


massive galaxy—its unusual l ocatio n far away in the outer halo 
(^proj ~ 54 knc; |[luxor et al.ll2014l) suggests that it formed in a 
dwarf galaxy like the LMC or Sgr. 

PA17’s low [Ca/Fe] and high [Mg/Ca] and [Ba/Eu] ratios indicate 
that PA 17 is more chemically similar to the LMC stars and 
clusters than to those associated with the MW (and presumably to 
those in M31), even when systematic offsets are considered. As 
discussed earlier, the [a/Fe] and [Ba/Eu] ratios indicate whether a 
cluster formed in an environment enriched by Type la supernova 
and AGB star products. Combined with its metallicity, PA17’s 
detailed chemical abundance ratios therefore indicate that it 
formed in a galaxy which had sufficient mass to have a fairly 
high star formation rate or to populate the high mass end of the 
IMF. However, its high [Mg/Fe] (and possibly its high [Na/Fe]) 
indicate that PA17 may have been enriched in ejecta from the 
most massive stars. As discussed in Sections 14.21 and 14.31 this 
suggests that PA17’s progenitor galaxy had more massive stars 
than a typical dwarf (i.e. it had a top-heavy IMF), was enriched 
by supernova ejecta from a rapidly rotating massive star, and/or 
was uniquely enriched from the products of a supernova due to 


inhomogeneous mixing in the host galaxy. 

Ultimately, PA17’s integrated abundances are chemically distinct 
from those of the metal-rich Galactic GC 47 Tuc, and are similar 
to Pal 1, the intermediate age LMC clusters, and the accreted Sgr 
clusters Pal 12 and Ter 7. 

H10 and H23: H10 and H23 are the third and second most metal- 
rich GCs in the sample of PAndAS GCs, respectively. Based on 
the same arguments as for PA 17, this indicates that the clus¬ 
ters likely originated in LMC- or Sgr-like dwarf galaxies. Un¬ 
like PA17, H23’s chemical abundance ratios are indistinguish¬ 
able from MW field stars and clusters and M31 GCs at the 
same [Fe/H], HlO’s [Ca/Fe] is mildly low while its [Mg/Ca] and 
[Eu/Ca] ratios are mildly high—this [Ca/Fe] ratio could indicate 
that H10 formed in a dwarf galaxy with a knee in the [Ca/Fe] 
vs. [Fe/H] relationship near [Fe/H] ~ —1.5, similar to LMC 
intermediate-aged clusters. 
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Figure 12. Syntheses of the 6645 A Eu II line. Grey regions show areas with uncertain HFS, while purple regions indicate uncertain molecular features. The 
red lines show the best fit, while the green and blue lines show ±lcr uncertainties. PA06 and PA56 have only upper limits for the Eu II abundance. 


5.1.2 PAndAS Clusters with [Fe/H] < —1.5 

It is more difficult to tag metal-poor GCs ([Fe/H] < —1.5) with 
chemical abundances, since the chemistries of metal-poor stars in 
dwarf galaxies and massive galaxies are not always significantly 
different. The analysis is further complicated by potential star-to- 
star Na, Mg, Ba, and Eu variations within the GCs. 

PA53 and PA56: These clusters are fairly close to each other 
in projection, and have very similar radial velocities and metal- 
licities; it is therefore possible that these clusters are physically 
associated with one another and were accreted from the same 
dwarf. Their discrepant Mg, Ba, and Eu ratios do not preclude the 
possibility that they are related, as the integrated abundances may 
be affected by star-to-star variations within the clusters. Both GCs 
are slightly enhanced in Na (with [Na/Fe]^ 0.4), which is likely 
a signature of the Na/O anticorrelation. PA56 also has a very low 
[Ba/Eu] (below the r-process only yields), similar to M15, which 
supports the idea that it may have star-to-star heavy element 
variations. Both clusters are metal-poor ([Fe/H] ~ —1-7) and 
have [Ca/Fe] ratios slightly lower than the MW field stars, similar 
to the Fornax field stars. 


PA54: Although PA54 is extremely close to PA53 in projection, 
their discrepant radial velocities imply that these clusters are not 
likely to be associated, nor is PA54 likely to be associated with 
PA56. However, PA54 has very similar abundance ratios as PA53 
and PA56, even with only upper limits on [Eu/Fe]. This indicates 
that PA54 may have formed in similar conditions. 

PA06: PA06 is the most metal-poor PAndAS GC in this analy¬ 
sis, though it is still more metal-rich than the Galactic GC M15. 
Again, PA06’s high [Na/Fe] indicates the presence of a Na/O 
anticorrelation—additionally, its low [Mg/Fe] hints at a Mg/Al 
anticorrelation. PA06 is similar to the metal-poor GCs in the MW 
(particularly Ml5) and in the dwarf galaxies: it is metal-poor and 
a-enhanced. The weakness of its spectral lines means that only an 
upper limit can be obtained for Eu II. 

5.2 Comparisons with M31 Halo Stars 

The presence of streams in the metal-poor density map of M31’s 
outer halo (see Figure [1} implies that some metal-poor (i.e. low 
mass) dwarf galaxies are currently being accreted. The chemical 
ratios of PA06, PA53, PA54, and PA56 are all consistent with an 
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Figure 13. Top: Comparisons of PAndAS clusters (red stars) t o MW field stars and dwarf galaxy field stars. Points and references are as in Figure[5a] The 
dashed red line shows the r-process-only limit from lBurris et al.l t200(] ). Bottom: Comparisons of PAndAS clusters (red stars) to MW field stars and GCs from 
the MW, M31, and various dwarf galaxies. Points and references are as in Figure [7bl 


accretion origin in at least one metal-poor dwarf galaxy. None of 
these GCs have been associated with streams based on their posi¬ 
tions. Given their lack of association with any streams, if PA06 and 
PA54 were accreted fro m dwarf satellites t hey may be remnants of 
ancient accretion events: flbata et al]j2QI4l ) estimate that ~ 58% of 
the metal-poor halo ([Fe/H] < —1.7) resides in a “smooth” com¬ 
ponent, which may include these GCs. 

Recall that PA53 and PA56 are located near the Eastern Cloud. 
Their similar kinematics and chemical compositions hint that PA53 
and PA56 are associated with each other—however, the numbers of 
stars around the clusters (Mackey et al. 2015, in prep.), their metal- 
licities, and their detailed chemical abundances indicate that PA53 
and PA56 were also likely associated with the dwarf galaxy that 
created the Eastern Cloud. In this case, the Eastern Cloud must turn 
back toward M31, encompassing PA53. This further suggests that 
the Eastern Cloud’s progenitor must have been massive enough to 
form GCs, though not sufficiently massive to produce a significant 
metal rich stellar population. 

The presence of the metal-rich Giant Stellar Stream (GSS), 
the associated H I gas, the outer halo stellar mass, and the num¬ 
ber of metal-rich GCs in the outer halo indicate that M31 likely 


experienced a minor merger with a Sgr or LMC-mass galaxy (e.g. 
[Fardal et ai1 l2013LlLewis et al.lEoTllBate et alj|2014l)1lbata etahl 
( 2014 ) find that 86% of the most metal-rich stars ([Fe/H] > —0.5) 
in the outer halo are associated with substructure from an accreted 
companion. Though this number decreases with metallicity, the 
fraction of accreted stars still remains high: 78% of the stars in 
PA17’s metallicity bin (—1.1 < [Fe/H] < —0.5) are expected to 
be associated with coherent, accreted structures, while 58% of the 
stars in H10 and H23’s metallicity bin (—1.7 < [Fe/H] < —1.1) 
are expected to be in streams. It is therefore likely that these com¬ 
paratively metal rich GCs were associated with dwarf satellites. 
H23 has been tent ative ly linked to a stream near the GSS (Stream 
D; see lVelianoski et al.l2014l) based on its position, though its kine¬ 
matics indicate that such an association is unlikely; neither H10 nor 
PA17 have been linked to any visible stellar streams, though HST 
imaging of fields around H10 and H23 reveal extremely metal-rich 
populations of field stars, suggesting that the GCs may be located 
on low surface brightness streams jRichardson et al.l2009l) . 


H10 is located near the SW Cloud, an overdensity of stars 
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to the southwest of M31 0 With photometry, iBate et all l20l4) 
estimate that the cloud hosts a metal-rich population ([Fe/H] ~ 
— 1.3), and that the progenitor galaxy was a fairly bright dwarf 
(My ~ —12, which is slightly fainter than For and Sgr). Three 
GCs (PA7, PA8, and PA14) appear to b e kinematically associated 
with each other and with the SW Cloud I Mac key et all2013Ll2014l : 
IVelianoski et alJl2014l : IBate et all 1 20141) . The SW Cloud does ex¬ 
tend to the sout heast, ancTHIO lies at the end of this extension (see 
Figure 2 in lBate et~aP[2014ll . HlO’s radial velocity (—352 km s -1 ) 
agrees with the stream’s radial velocity and that of the GCs, and 
even agrees with the velocity gradient noted by Bate et al. (where 
the northernmost GC, PA7, is moving toward the MW faster than 
PA8 and PA17). Furthermore, HlO’s [Fe/H] ~ —1.4 agrees very 
well with the Cloud and its GCs, and its detailed abundances sup¬ 
port the idea that it could have originated in a dwarf galaxy with 


16 Note that though PA17’s proximity to PA14 hints at an association with 
the SW Cloud, its dra stically different radial velocity (—260 km s 1 ) from 
PAM (—363 km s —1 VeljanosHet al ill) 111 ) and high [Fe/H] (—0.93 com¬ 
pared to —1.30: lMackev et~aI]|2oT3i) ~make it unlikely that PA17 was asso¬ 
ciated with the SW Cloud progenitor. 


the mass of For or Sgr. Thus, it is likely that H10 is associated with 
PA7, PA8, PAM, and the SW Cloud. 


As mentioned above, H23 has been spatially associated with 
Stream D based on its position, though its velocity disagrees with 
the other GCs in the stream. Its metallicity, [Fe/H] = —1.1, now 
further suggests that H23 is not likely to be associated with that 
stream, since Stream D is primar ily com posed of metal-poor stars 
with [Fe/H] < —1.1 dlbata et alj|2014l) . It is more likely that H23 
came f rom the progenitor of the GSS, along with PA17. lFardal et al] 
d2013l) perform N-body simulations to reproduce the GSS and other 
stellar debris. These simulations indicate that a massive progenitor 
was accreted during multiple orbits, culminating in a recent final 
accretion (~ 760 Myr ago). Many of the metal-rich PAndAS GCs 
might have been stripped early on, and may no longer be associated 
with any bright, coherent substructure. Given their metallicities and 
other chemical abundances, it is possible that H23 and PA 17 both 
formed in the GSS progenitor. If H23 and PA 17 were born in the 
same dwarf galaxy, their [Ca/Fe] ratios clearly imply a “knee” at 
[Fe/H] ~ —1.1 to -1.3, similar to the LMC stars and clusters. 
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6 CONCLUSIONS 

This paper has presented integrated Fe, Na, Mg, Ca, Ti, Ni, Ba, and 
Eu abundances of seven outer halo M31 GCs, five of which were 
discovered in the Pan-Andromeda Archaeological Survey. These 
detailed IL chemical analyses of PAndAS clusters have identified 
GCs in an extragalactic system that may have been accreted from 
dwarf galaxies. Detailed investigations of the chemical abundance 
ratios of individual targets such as this are only possible with high 
resolution spectroscopy. 

The PAndAS cluster abundances suggest that these outer halo 
M31 GCs may have been accreted from multiple dwarf galaxies: 

• The metal-rich GC PA 17 is chemically more similar to the 
LMC stars and clusters than to MW field stars and clusters, 
suggesting that it originated in an LMC-like dwarf galaxy. 
H23’s abundances are indistinguishable from MW field stars 
and clusters, though its metallicity and location suggest it may 
have originated in a massive dwarf galaxy. PA 17 and H23 may 
have been accreted along with the progenitor of the metal-rich 
Giant Stellar Stream. 

• The intermediate metallicity GC H10 has a location, metallicity, 
and radial velocity that agree well with the SW Cloud and its 
GCs, which may have originated in a Sagittarius or Fornax sized 
progenitor l lBate et al.1120141) . HlO’s chemical abundance ratios 
support the suggestion that it may have formed in a dwarf galaxy. 

• PA53, PA54, and PA56 have abundances and metallicities that 
are more typical of an intermediate mass dwarf galaxy like 
For. This suggests that they are currently being accreted from 
at least one metal-poor dwarf galaxy and could be associated 
with one or more of the coherent, metal-poor streams observed 
in PAndAS. Based on their chemistries and radial velocities 
PA53 and PA56 could be associated and are likely to have been 
accreted from the dwarf galaxy that created the Eastern Cloud. 
Despite its proximity in projection, it is unlikely that PA54 is 
associated with either PA53 or PA56. 

• PA06’s metallicity makes it a difficult target for chemical 
tagging analyses, since the chemistries between dwarf and 
massive galaxies have likely not had sufficient time to diverge 
at [Fe/H] ~ —2. Several of its integrated abundances do not 
agree with MW or dwarf galaxy stars, suggesting that strong 
star-to-star chemical variations are present in the cluster. 

Thus, this chemical tagging analysis is consistent with the observa¬ 
tion that M31’s outer halo GC system is being built up by accretion 
of dwarf satellites. In addition, this detailed abundance analysis 
provides additional information on the nature of these progenitor 
systems that is, at present, very difficult to obtain from the field 
populations. 
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APPENDIX A: RGB SLOPE VS. METALLICITY: A 
CALIBRATION WITH GALACTIC GLOBULAR CLUSTER 
OPTICAL CMDS 


It is well established that the slo pe of a cl uster ’s RGB is correlated 
with cluster metallicity (e.g. ISaraiedini|[l994l) . This appendix dis¬ 
cusses the calibration of RGB slope vs. [Fe/H], utilizing high qual¬ 
ity HST CMDs of Galactic GCs from the ACS Survey of Galactic 
Globular Clusters iSaraiedini et ai]l2007l: lAnderson et al.ll2008l) FI 
For this calibration 15 Galactic GCs with high quality CMDs and 
well-populated RGBs that span a wide range in [Fe/H] were se¬ 
lected; priority was given to clusters with low foreground redden¬ 
ing. Using the origin al F606W and F814W H ST magnitudes, the 
distance moduli from lVandenBerg et alj d2013l) were used to over¬ 
plot clusters, and distance moduli and reddening values were ad¬ 
justed so that cluster HBs overlapped. The average RGB colours 
were then determined at two magnitudes (My = 0 and Mv = 
—2), and RGB slopes were calculated. These points are shown in 
Figure [AT] while the slopes are listed in Table [AT] 

Uncertainties in RGB slope were estimated based on the un¬ 
certainties in colour at the two magnitudes. The metal-poor clus¬ 
ters have small uncertainties in RGB colour, though this translates 
into large uncertainties in slope because the RGBs are steeper. The 
metal-rich clusters have very uncertain RGB colours, but because 
the RGBs are flatter the uncertainty in slope is much lower. Note 
that at higher [Fe/H] (> —0.7) a linear fit to the RGB is no longer 
a good approximation to the actual shape of the RGB—thus, this 

calibration is likely to break down at the metal rich end._ _ 

For each cluster the [Fe/H] ratios from lCarretta et all (1200% ) 
were adopted; [Fe/H] uncertainties of 0.05 dex were assumed. Fig¬ 
ure IA2I shows the relationship between RGB slope and cluster 
[Fe/H], The linear least squares fit to the points is shown with the 
solid line, while the dashed lines show the uncertainty in the fit. 
The trend is clear: metal-poor GCs have steeper RGBs than more 
metal-rich GCs. Non-standard chemical abundance mixtures (in, 
e.g., C, N, O or a-elements like Mg or Si) can affect the s hape of 
the RGB (e.g. ISalaris etliill 1 9931 : IVandenBerg et al.|[2012l) . which 
would affect where a cluster falls in the plot. The adopted [a/Fe] 
ratio is particularly important: the qualitative effects of [a/Fe] on 
the F606W, F814W RGB slope are illustrated in Figure I A3I u sing 
the Victoria-Regina isochrones from VandenBerg_etal. (2014^_and 
the colour transformations of ICasagrande & VandenBergl d2014h . 
Figure [3al shows that at a given [Fe/H], clusters with [a/Fe] m 0 
have steeper RGBs than clusters with [a/Fe] = +0.4; this offset 
is greater at higher metallicity. Furthermore, Figure [3bl illustrates 
that the slope difference can be mimicked by lowering [Fe/H] by 
0.3 dex while maintaining [a/Fe] = +0.4 (note that in Figure [3bl 
the dotted isochrones have been shifted by 0.02 magnitudes so that 
both the a-enhanced and a-normal isochrones can be seen). The 
blue dot-dashed line in Figure IA2I shows the effects of lowering 
[a/Fe] by 0.4 dex, which is approxim ated by lowering the [Fe/H] 
by 0.3 dex (also see lSalaris et al JlT993h . 

The RGBs of the PAndAS clusters are shown in Figure lA4l 
while the RGB slopes and spectroscopic metallicities are shown in 
FigurelTbl The photometric uncertainties, sparseness of the CMDs, 
and difficulties in determining distance moduli (particularly for 
PA56) make the RGB slopes more uncertain for the PAndAS clus¬ 
ters, compared to the Galactic GCs. However, all PAndAS GCs 
agree with the MW RGB slope relation in Figure lA2l within their 


Table Al. Galactic GC RGB slopes. 


Cluster 

[Fe/H] 

RGB slope 

M92 

-2.35 + 0.05 

-11.429+1.172 

M15 

-2.33 + 0.05 

-11.765+1.238 

M53 

-1.93 + 0.05 

-10.929+1.077 

NGC 5286 

-1.70 + 0.05 

-10.870+1.066 

M2 

-1.66 + 0.05 

-9.346+0.799 

M13 

-1.58 + 0.05 

-9.217+0.778 

M3 

-1.50 + 0.05 

-9.091+0.758 

M5 

-1.33 + 0.05 

-8.299+1.181 

NGC 1261 

-1.27 + 0.05 

-7.937+1.087 

NGC 1851 

-1.18 + 0.05 

-7.491+0.976 

NGC 6362 

-1.07 + 0.05 

-6.557+0.760 

47 Tuc 

-0.76 + 0.05 

-3.766+1.031 

NGC 6652 

-0.76 + 0.05 

-4.090+1.187 

M69 

-0.59 + 0.05 

-3.883+1.086 

NGC 5927 

-0.29 + 0.05 

-7.491+0.595 


References: Cluster [Fe/H] ratios are from Carr e tta et al] ll2009b ). 
Table A2. PAndAS GC RGB slopes. 


Cluster 

[Fe/H] 

RGB slope 

PA06 

-2.06 + 0.10 

-12.0482+2.658 

PA54 

-1.84 + 0.10 

-11.173+1.405 

PA56 

-1.73 ± 0.10 

-10.811+2.982 

PA53 

-1.64 + 0.10 

-10.526+1.238 

H10 

-1.40 + 0.10 

-8.621+0.813 

H23 

-1.12 ± 0.10 

-7.353+1.268 


References: Cluster [Fe/H] ratios are from the [Fe I/H] abundances in 
Table [6] 

errors. This relationship illustrates that the spectroscopic metallici¬ 
ties are consistent with the observed RGB slopes in the HST CMDs. 
All of the clusters in Figure [4b] lie above the MW relation closer to 
the low [a/Fe] value, hinting at the possibility of low [a/Fe] ratios 
in the PAndAS targets. 


17 http://www.astro.ufl.edu/~ata/public_hstgc/ 
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Figure Al. F606W, F814W photometry of Galactic GC RGBs, from the ACS Globular Cluster Survey fcaraiedini et al.|[2007l ; lAnderson et alj|2008l) . The 
cluster HBs were aligned, and the colours at two magnitudes ( My = 0 and My = —2, shown as the horizontal dashed lines) were d etermined for each 
cluster; these values are shown as red crosses. The clusters are offset in the plot, and are ordered by metallicity (from lCarretta et all2009al) . 



RGB slope 


Figure A2. RGB slope versus [Fe/H] for the Galactic GCs. Slopes are determined using the points in Figure IaTI slope errors were estimated based on the 
uncertainties in these RGB colours. 
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Figure A3. Victoria-Regina F606W, F814W isochrones from lVandenBerg et al.H2014l) . Left: Metal poor ([Fe/H] = —1.7, in blue) and metal rich ([Fe/H] = 
— 1.0, in red) isochrones at two different values of [a/Fe]. The solid lines show [a/Fe] = +0.4, while the dashed lines show [a/Fe] = 0.0. Isochrones 
with low [a/Fe] have steeper RGB slopes at a fixed metallicity. Right: Isochrones illustrating that the steeper slopes from [a/Fe] = 0.0 isochrones can be 
reproduced by lowering [Fe/H] by 0.3 dex while maintaining [a/Fe] = +0.4. The dashed lines show the same isochrones as in Figure[3a] while the dotted 
lines show isochrones with [a/Fe] = +0.4 and A[Fe/H] = —0.3 dex. The dotted isochrones are offset by +0.02 magnitudes so that both isochrones can be 
seen. This shows that an [Fe/H] offset of -0.3 dex can be used in the RGB slope calibration to approximate the effects of lower [a/Fe]. 



Figure A4. F606W, F814W HST photometry of PAndAS GCs, from Mackey et al. (in prep.). The cluster HBs were aligned, and the colours at two magnitudes 
(My = 0 and My = — 2, shown as the horizontal dashed lines) were determined for each cluster; these values are shown as red crosses. The clusters are 
offset in the plot, and are ordered by metallicity. 
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